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ABSTRACT 



A method is presented to derive electron temperatures and densities of plan- 
etary nebulae (PNe) simultaneously, using the observed hydrogen recombination 
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spectrum, which includes continuum and hne emission. By matching theoreti- 
cal spectra to observed spectra around the Balmer jump at about 3646 A, we 
determine electron temperatures and densities for 48 Galactic PNe. The elec- 
tron temperatures based on this method [hereafter Te(Bal)] are found to be 
systematically lower than those derived from [O III] A4959/A4363 and [O III] 
(88/im + 52/im)/A4959 ratios [hereafter, Te([0 III]„a) and Te([0 III]f„)]. And the 
electron densities based on this method are found to be systematically higher than 
those derived from [O II] A3729/A3726, [S II] A6731/A6716, [CI III] A5537/A5517, 
[Ar IV] A4740/A4711 and [O III] 88//m/52//m ratios. These results suggest that 
temperature and density fluctuations are generally present within nebulae. The 
comparison of ^^([O IIl]na) and T'e(Bal) suggests that the fractional mean-square 
temperature variation (t^) has a representative value of 0.031. A majority of 
temperatures derived from the [O III] (88;um -|- 52fim) / A4959 ratio are found to 
be higher than Tc([0 IIl]na), which is attributed to the existence of dense clumps 
in nebulae - those [O III] IR fine-structure lines are suppressed by collisional dc- 
excitation in the clumps. By comparing Te([0 III]fn), Te{[0 IIl]na) and Te(Bal) 
and assuming a simple two-density-component model, we find that the filling 
factor of dense clumps has a representative value of 7xl0~^. The discrepancies 
between Te{[0 IIl]na) and Te(Bal) arc found to be anti-corrclatcd with electron 
densities derived from various density indicators; high-density nebulae have the 
smallest temperature discrepancies. This suggests that temperature discrepancy 
is related to nebular evolution. In addition, He/H abundances of PNe arc found 
to be positively correlated with the difference between Te([0 III]na) and Te(Bal), 
suggesting that He/H abundances might have been overestimated generally be- 
cause of the possible existence of H-deficient knots. Electron temperatures and 
densities deduced from spectra around the Paschen jump regions at 8250 A are 
also obtained for four PNe, NGC 7027, NGC 6153, M 1-42 and NGG 7009. Elec- 
tron densities derived from spectra around the Paschen jump regions are in good 
agreement with the corresponding values derived from spectra around the Balmer 
jump, whereas temperatures deduced from the spectra around the Paschen jump 
are found to be lower than the corresponding values derived from spectra around 
the Balmer jump for all the four cases. The reason remains unclear. 

Subject headings: atomic processes - planetary nebulae: general - ISM: abun- 
dances 
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1. Introduction 

Accurate determinations of metal abundances in planetary nebulae (PNe) are important 
for understand the chemical evolution of galaxies and the nucleosynthesis and mixing pro- 
cesses in low- and intermediate-mass stars. A long-standing problem in nebular abundance 
studies has been that the heavy-element abundances derived from optical recombination 
lines (ORLs) are systematically higher than those derived from coUisionally excited hnes 
(CELs; see Liu 2001, 2003 for recent reviews). In extreme cases, the discrepancies exceed a 
factor of 10. Given that GEL abundances are sensitive to electron temperature and density, 
knowledge of temperature and density variations is an essential ingredient in understanding 
the abundance discrepancy problem. Information on temperature and density variations 
can be obtained by comparing results derived from plasma diagnostics which behave dif- 
ferently as temperature and density vary. Traditionally, nebular electron temperatures and 
densities are derived from GEL ratios. However, for PNe including high-density regions 
(A^e> lO^cm-^), such as IG 4997 (Hyung et al. 1994), Mz 3 (Zhang & Liu 2002) and M 2-24 
(Zhang & Liu 2003), GELs are strongly suppressed by coUisional de-excitation, and so de- 
pending on their critical densities, they may become unusable in probing such high-density 
regions. On the other hand, hydrogen recombination lines have very high critical densities 
(> lO^cm"^). Consequently, the hydrogen recombination spectrum provides a powerful tool 
to probe high-density regions. Comparison of electron temperatures and densities derived 
from the hydrogen recombination spectrum (lines plus continuum) with those derived from 
CELs should enable us to quantify temperature and density variations in nebulae. 

By measuring the Balmer discontinuity, Liu & Danziger (1993) determined electron 
temperatures in 14 PNe and found that they are systematically lower than those derived from 
the [O III] nebular-to-auroral forbidden line ratio. This can be explained by the presence of 
large temperature fluctuations since the [O III] A4363 auroral line of higher excitation energy 
tends to be emitted in higher temperature regions than the [O III] AA4959, 5007 nebular hnes 
of lower excitation energy, as suggested by Peimbert (1967). Recent reviews on temperature 
fluctuations of gaseous nebulae are given by Torres-Peimbert (2003) and references therein. 
However, no photoionization models yield such large temperature fluctuations suggested by 
observations. Moreover, temperature fluctuations are unable to explain the large disparity 
between the ORL and CEL abundances observed in some PNe (Liu 2001). Instead, it is 
found that the hypothesis that PNe may contain clumps of different temperature, density 
and chemical composition, embedded in diffuse 'normal' material, can explain many of the 
observed patterns (Liu et al. 2000). Peimbert's concept of temperature fluctuations alone 
may not be a realistic description of the physical conditions in a nebula, although it has 
proven to be a convenient and useful way to characterize the complexities of real nebulae. 
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It is expected that there are also variations of electron density inside a PN. Rubin (1989) 
pointed out that electron densities derived from various plasma diagnostics should correlate 
positively with the critical densities of the diagnostic lines involved if there are large density 
variations within a nebula. Liu ct al. (2001a) found that electron densities derived from the 
[O III] 88/im/52/im ratios arc generally lower than those derived from the optical [Ar IV] 
and [CI III] doublet ratios, suggesting the presence of high density regions in PNe, where the 
52- and 88-/im lines are suppressed by coUisional de-excitation due to their relatively low 
critical densities. Viegas & Clegg (1994) showed that the discrepancies between Te{[0 llljna) 
and Te(Bal) can be explained by the presence of condensations in nebulae and gave an 
analytical method of estimating the filling factor of condensations. Our detailed studies of 
the peculiar PNe Mz 3 and M 2-24 suggest that the presence of high-density regions in some 
PNe can significantly affect abundance determinations (Zhang & Liu 2002, 2003). However, 
density inhomogcncities alone cannot explain the large disparity between the ORL and CEL 
abundances either, unless the condensations are also H-deficient (Liu et al. 2000). 

Liu et al. (2000) found that the discrepancy between the ORL and CEL abundances 
correlates with the difference between the temperatures derived from the [O III] collisional 
lines and from the Balmer jump of hydrogen recombination spectrum, thus lending strong 
support to the argument that some extremely cold H-deficient knots may co-exist with normal 
nebular material in many PNe. However, it is difficult to explain the origin and evolution of 
such H-deficient knots with the current theories of element production in intermediate-mass 
stars. 

In this paper, we present determinations of electron temperatures and densities for a 
large sample of PNe based on the analysis of nebular hydrogen recombination spectrum. Sec- 
tion 2 describes the sample. In Section 3, we describe a new approach to deriving electron 
temperature and density simultaneously from the observed nebular hydrogen recombination 
spectrum (lines plus continuum) and present the results. In Section 4, we compare results 
derived from a variety of temperature and density indicators and discuss the significance of 
variations in electron temperature and density in PNe. A discussion of He/H abundance 
determinations is also presented in this Section. In Section 5, we present electron temper- 
atures and densities derived from the Paschen discontinuity and decrements. A summary 
then follows in Section 6. 
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2. The sample 

2.1. Observations and data reduction 

48 Galactic PNe (37 disc PNe and 11 bulge PNe) are analyzed in the current work. The 
spectra were obtained in a series of observing runs using telescopes at ESO La Silla and La 
Palma. 

14 objects were observed using the ISIS double long-slit spectrograph mounted on the 
4.2 m William Herschel Telescope (WHT) at La Palma Observatory in 1996 July-August and 
in 1997 August. A Tek 1024 x 1024 24//m x 24//m charge coupled device (CCD) was used, 
which yielded a spatial sampling of 0.3576 arcsec per pixel projected on the sky. In 1996, a 
600gmm~^ and a 316gmm~^ grating were used for the Blue and Red Arm, respectively. In 
1997, they were replaced by a 1200 gmm~^ grating. 

23 objects were observed with the ESO 1.52 m telescope using the Boiler & Chivens 
(B&C) long-slit spectrograph. The spectra were obtained during runs in 1995 July, 1996 July, 
1997 February and 2001 June. In 1995 the detector was a Ford 2048 x 2048 15fim x 15fim chip, 
which was superseded in 1996 by a thinned UV-enhanced Loral 2048 x 2048 15;um x 15;um 
chip of much improved quantum efficiency. The B&C spectrograph has a useful slit length 
of about 3.5arcmin. In order to reduce the CCD read-out noise, the CCD was binned by 
a factor of two along the slit direction, yielding a spatial sampling of 1.63 arcsec per pixel 
projected on the sky. A slit width of 2 arcsec was used throughout the ESO 1.52 m runs. 

11 objects were observed using the Intermediate Dispersion Spectrograph (IDS) on the 
2.5m Isaac Newton Telescope (INT) at the Observatorio del Roque de los Muchachos, on La 
Palma, Spain in 2001 August. A slit 4 arcmin long and 1 arcsec wide was used for nebular 
observations. The R1200B grating was used to cover wavelengths from 3500-5000 A at a 
spectral resolution of 0.95 A, determined from the FWHM oi the cahbration arc lines, while 
the R300V grating, together with a GG385 order sorting filter, was used to cover wavelengths 
fi-om 3800-7700 A at a resolution of 4.0 A FWHM. 

All the spectra were reduced using MIDAS^ or IRAF^ following the standard pro- 
cedure. The spectra were bias-subtracted, fiat-fielded and cosmic-rays removed, and then 

wavelength calibrated using exposures of a calibration arc lamp. The spectra were flux- 
calibrated by observing spcctrophotometric standard stars. 1-D spectra were obtained by 
integrating along the slit. Using the Galactic reddening law of Howarth (1983), the log- 



-'^ MIDAS is developed and distributed by the European Southern Observatory. 
^IRAF is developed and distributed by the National Optical astronomy Observatory. 
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arithmic extinction at H/3, c = log/(H/9)/F(H/3), was derived from the observed Balmer 
decrements, the Ha/H/?, H7/H/3 and H5/H/3 ratios. A mean value for each nebula was used 
to deredden the spectra. The object names and the extinction constants are given in Table 1. 

Amongst the PNe in our sample, far-infrared spectra for 21 of them have been obtained 
using the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory 

{ISO). Fluxes of far-IR fine structure lines derived from these spectrum have been published 
by Liu et al. (2001a) except for three objects, Hu 1-2, NGC 6210 and NGC 6818. For the 
latter three PNe, we measured far-IR line fluxes using spectra retrieved from the ISO data 
archive. In order to normalize the ISO line fluxes to H/?, the absolute H/3 fluxes tabulated 
in Cahn et al. (1992) were used. 



2.2. He+/H+ and He2+/H+ abundance ratios 

In order to calculate theoretical nebular spectra, accurate determinations of the He+/H+ 
and He^+/H+ abundance ratios are necessary (see Section 3). Helium ionic abundances were 
derived from He I and He II recombination lines, using 

Ar(He'+) _ /a A anp . . 

N{m) Inp Xnp «He^+ ' ^ ^ 

which have only a weak dependence on the adopted electron temperature and density. 

He^/H"*" abundance ratios were calculated from the strong line He IA5876 under the 
assumption of Case A. The effective recombination coefficient 0^0+7. derived from the 
formulae presented by Benjamin et al. (1999). He^+/H+ abundance ratios were calculated 
from the He IIA4686 line, using the effective recombination coefficient 01-^^2+ of Storey & 

4686 

Hummer (1995). The results are presented in Table 1. 



2.3. Temperature and Density determined from far-IR and optical forbidden 

line ratios 

For the purpose of comparison, various forbidden line ratios are also used to determine 
electron temperatures and densities. For densities, results derived from the following diag- 
nostic ratios were used, [O III] 52 //m/88 //m, [O II] A3729/A3726, [S II] A6716/A6731, [CI III] 
A5517/A5537 and [Ar IV] A4711/A4740. Electron densities from the [O III] 52-//m/88-//m 
ratio are taken from Liu et al. (2001a) except for NGC 6210 and NGC 6818, for which we 
obtain the values of logA^e — 3.19 and 2.86 (cm~^), respectively, using our own measured 
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line fluxes. Densities deduced from the optical [O II] A3729/A3726, [S II] A6716/A6731, 
[CI III] A5517/A5537 and [Ar IV] A4711/A4740 doublet ratios are taken from Wang et al. (in 
preparation). A detailed discussion of these density indicators is given in their paper. 

Electron temperatures derived from the [O III] (88//m + 52//m)/A4959 and [O III] 
A4959/A4363 ratios together with their uncertainties are presented in Table 1. For the 
former, values were derived assuming electron densities derived from the [O III] 52-/im/88- 
fim line ratio. For the latter, the average densities derived from the [CI III] A5517/A5537 
and [Ar IV] A4711/A4740 doublet line ratios were adopted, given that C1+ and Ar^+ have 
ionization potentials of 23.8 cV and 40.7 cV, respectively, comparable to the value of 35.1 eV 
for O"*", and therefore anticipating that the [O III] lines arise from similar ionized regions 
as [CI III] and [Ar IV] lines. As the density increases, the resultant temperatures become 
increasingly sensitive to the adopted density. When the density reaches a value higher than 
10^ cm-^ the [O III] (88/im + 52//m)/A4959 and [O III] A4959/A4363 ratios are no longer 
useful temperature diagnostics. For example, high density regions have been detected in 
M 2-24 and IC 4997 (Zhang & Liu 2003; Hyung et al. 1994). For both PNe, temperatures 
deduced from the [O III] A4959/A4363 ratio have large uncertainties. 

3. Electron temperature and density determined from the hydrogen 

recombination spectrum 

The H I recombination spectrum consists of continuum and line emission. The continu- 
ous spectrum has a jump at the limit of the Balmer series at 3646 A. High-order Balmer lines 
(n — 2, n > 25) converge at the Balmer discontinuity and merge with the continuum jump 
at shorter wavelength. The size of the Balmer jump is a function of electron temperature. 
On the other hand, intensities of high-order Balmer lines (n — > 2, n > 10) relative to H/3 
are sensitive to electron density. Therefore, the hydrogen recombination spectrum near the 
Balmer jump region provides valuable temperature and density diagnostics. Our method 
and results of the electron temperature and density determinations based on this technique 
are presented below. 

3.1. Theoretical spectra 

Several processes contribute to the continuum emission from a nebula: free-free and free- 
bound emission and the two-photon decay. Apart from hydrogen, the second most abundant 
element, helium also contributes significantly to the first two processes. In addition, emission 
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from the central star and local dust can in principle contaminate the observed continuum 
spectrum. In the following, a detailed calculation including all the processes is presented. 

The free-free (or bremsstrahlung) emission, which is emitted when an free electron 
interacts but does not combine with a positively charged ion such as H+, He+ and He^+, 
dominates the nebular continuum in the infrared and radio region. The emission coefficient 
of this process is given by Brown & Mathews (1970), 

where iVj and are ion and electron densities respectively, Z is the nuclear charge, Tg is 
the electron temperature, hfQ denotes the Rydberg energy and Qff is the free-free Gaunt 
factor. Here gff is computed using the method outlined by Hummer (1988). 

In the optical region, free-bound emission is the main contributor to the continuum. 
According to Brown & Mathews (1970), the emission coefficient of continuous emission due 
to recombinations of H+, He"^ or He^"*" with electrons is 

Air Trmekle 

n=ni 1=0 

where huni is the ionization potential of the {n,l) state, rii is the lowest state that can 
contribute to the free-bound emission at the given frequency u and a„i(i/) is the photoion- 
ization cross-section for the state, which is computed using the method described by Storey 
& Hummer (1991). 

Another important source of continuous emission is two-photon decay of the 2^5" level 
of the hydrogen atom (two-photon emission from He+ and He° is neghgible for nebulae of 
cosmic compositions; Brown & Mathews, 1970). Single photon transition from the level 
to the ground I'^S is forbidden. Radiative decay of the level to the ground can only 
occur via two photon emission, in which the electron jumps to a virtual level, which can lie 
in anywhere between n = 1 and n = 2 levels, and then to the ground state, emitting two 
photons in the process. As a result, continuum emission longward of Lya will occur. The 
emission coefficient of the two-photon process is given by 

ju,2q = ■^n22shyA{y) [erg cm~^ s~^ Hz~^], (4) 

where y = u/ui^ya and A{y)dy is the probability per second for emitting a photon in the 
interval dy. Approximately, the function A{y) is given by the following expression 

Aiy) = 202.0{y(l - y){l - [4y(l - y)]°«} + 0.88[y(l - y)]^-^^[4y(l - y)]°«}s-^ (5) 



-9- 



(Kwok 2000). The transition probabihty for the two-photon decay is A225 125 = A{y)dy — 
8.23 s~^. The equihbrium population of the 2^5" state is given by 

N N nH^ 

^2^5 = 1 , AT P , AT e (6) 

^225,125 + ^^p522 5,22p + ^^e922s,22p 

(Osterbrock 1989), where Np is the proton density, CK^lg is the effective recombination coeffi- 
cient for populating 2^5*, and (l\2g2'^p and q'225 22p ^^e collisional de-excitation transition 
rates of 2^5* to 2^P states by electrons and protons, respectively. Note that the collisional 
transition from 2^P to 2^5* has been ignored in equation (6). Under Case A conditions, this 
is a good approximation. Under Case B conditions, however, the population of the 2^P state 
is comparable to that of the metastable 2^5" level since Lyct is trapped between the 2^P and 
I'^S states. Consequently, there is considerable colhsional excitation from 2^P to 2^5", which 
enhances the two-photon emission. Therefore, the emission coefficient of two-photon process 
given by equation (4) can be significantly underestimated. However, this hardly affects our 
resultant electron temperatures and densities since an additional continuum (see below), is 
added to fit our observed spectrum in order to simulate direct and scattered stellar light. 
For a limited wavelength region, such an approach is sufficient to account for any excess 
emission from the two-photon process. 

Finally, direct or scattered light from the central star can contaminate the observed 
continuum level. For our observations, the slit position was normally chosen to sample the 
brightest parts of the nebulae, sometimes passing through the central star. Some spectra were 
derived by scanning across the nebular surface. It is quite difficult to evaluate the level and 
spectral energy distribution of the contaminating stellar continuum. As central stars of PNe 
are relatively hot (>20,000 K), in our calculation, a Rayleigh- Jeans approximation for a Black 
Body (ja,rj ^~'^) is adopted to simulate the stellar continuum. This stellar component also 
partially compensates for the underestimated two-photon continuum, as mentioned above. 
For simplicity, we shall call this term 'stellar continuum' hereafter. Continuum emission 
from dust grains peaks in the infrared and is negligible in the optical region. Thus it has 
been neglected in our calculation. 

Apart from continuum emission, recombination line emission from bound-bound tran- 
sitions of H°, He° and He"*" is also included. The emission coefficient of a recombination line 
is given by 

3k,bb = ^NeNia'/^huk [erg cm'^ s"^], (7) 

where a^^'^ is the effective recombination coefficient of the line. For H I and He II recombi- 
nation lines, emission coefficients are determined as in Hummer & Storey (1987). Here the 
considered maximum principle quantum number of is n = 500. For He I recombination lines. 
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emission coefficients are taken from Benjamin et al. (1999) and Brocklefiurst (1972). All line 
profiles are assumed Gaussian with a FWHM 

= ^Jalt + al A, (8) 
where (Tk,t is the FWHM caused by thermal broadening, given by 

a,, = 2^X,{^f/' A, (9) 

where m is the mass of the atom, and (Jq is the FWHM caused by other broadening effects, 
such as nebular expansion, turbulence, seeing conditions and instrument. For a typical 
PN observed here, ao ^ Sck^t- Line emission coefficients are thus converted into emission 
coefficients per unit wavelength by 

jx,bb = Y] — ^exp(— 4^)jfe,66 [erg cm'^ g-^ A"^]. (10) 

Just longward of the Balmer jump, high-order Balmer lines converge and merge into a 
continuum-like feature. The slope of the feature is sensitive to electron density, and thus 
provide an excellent density diagnostic. The region is contaminated by the Ne II A3664 line 
of Multiplet V 1. However, the weak Ne II recombination line has a relatively narrow line 
width compared to the profile of the continuum-like feature and consequently, cannot cause 
any difficulty in our analysis. 

Summing up contributions from all the processes discussed above, the flux at a given 
wavelength is given by 

jx,tot — 3\,2q + 3\,RJ + JA,6/(X) jA,66(X)] [erg cm ^s ^A ]. (11) 

X=HO,He°,He+ 

Equation (11) shows that, in order to synthesize the nebular spectrum near the Balmer 
jump wavelength region, the required input parameters include nebular electron temperature 
and density, the magnitude of contaminating stellar continuum, He+/H+ and He^+/H+ ionic 
abundance ratios and cTq. In general, decreasing (increasing) of the nebular electron tem- 
perature will steepen (flatten) the general shape of the theoretical continuum and increase 
(decrease) the Balmer discontinuities; decreasing (increasing) of electron density will steepen 
(flatten) the continuous feature redwards of 3646 A and decrease (increase) the intensities of 
high-order Balmer lines relative to a low-order line, such as H/3, as shown in Fig. 1. Accord- 
ingly, matching the dereddened observed spectrum to those calculated from equation (11) 
by varying the electron temperature and density, we can determine the electron tempera- 
ture and density simultaneously provided other parameters, such as He+/H+ and He^+/H+ 
abundance ratios and Uo are known. 
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3.2. Results 

Electron temperatures and densities for 48 PNe are derived by fitting the synthetic to the 
observed spectra. For synthesis, the required He+/H+ and He^+/H+ abundance ratios have 
been observationally determined, and Go was derived by fitting the hne profiles of strong 
H I recombination lines deconvolved with the thermal broadening, as listed in Table 1. 
Three free parameters, nebular electron temperature and density and the fraction of the 
scattered Ray leigh- Jeans stellar continuum, are optimized simultaneously, using least Chi 
square method to match the observed spectra. For the optimization, the fitted data points 
were specified explicitly by eye, within the wavelength range from 3550-3850 A, but excluding 
regions of strong (forbidden) lines. The results are partially shown in Fig. 2 together with 
the residuals of the fits. 

In Table 2, we report the derived electron temperatures and densities along with their 
estimated errors. The Balmer jump and decrement of the H I recombination spectrum have 
previously been used to determine nebular electron temperatures and densities in four PNe, 
M 2-24 (Zhang & Liu 2003), NGC 6153 (Liu et al. 2000), M 1-42 and M 2-36 (Liu et al. 
2001b). For the four PNe, the results derived here are in good agreement with the earlier 
results within the errors. 



3.3. Error analysis 

The errors caused by uncertainties in the input parameters He^/H+ and He^+/H^ abun- 
dance ratios are negligible due to the small measurement errors of strong He I and He II lines 
and their small contribution to the continuum emission compared to hydrogen (~ 10%). Er- 
rors in (To can also be ignored as this quantity is well determined by fitting strong hydrogen 
hnes. 

The complicated velocity fields within a nebula can cause different lines to have different 
deviations from their laboratory wavelengths, and thus can affect the quality of line-fitting 
made under the simple assumption of a uniform velocity for all ions and lines. Line- fitting can 
also be affected by blends of unknown lines and deviations of line profiles from Gaussian. 
Inspection of Fig. 2 shows that the observed intensities of the He I 2s^S-5p^P° A3614 
lines are generally lower than its predicted values. Liu et al. (2001b) attributed this to 
the destruction of He I Lyman photons by photoionization of H° or by absorption of dust 
grains. The wavelength region near the Balmer jump also contains a few strong O III Bowen 
fluorescence lines which are not considered in our simulations, thus giving large residuals 
at the corresponding wavelengths. However, these effects do not contribute to errors in the 
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derived electron temperatures and densities since these wavelengths have been excluded in 
the fitting. 

The observed spectra plotted in Fig. 2, were de-reddened using extinction constants 
listed in Table 1. All spectra were normalized to H 11 A3770 in order to minimize errors in 
derived temperatures and densities caused by uncertainties in the reddening correction and 
flux calibration, and taking the advantage of the small wavelength difference between the 
Balmer jump and H 11. 

The stellar continuum may have a spectral energy distribution that deviates from the 
Ray leigh- Jeans approximation adopted here. However, this continuum is expected to be 
smooth and vary only slightly over such a short wavelength range of fitting. Therefore, we 
have also ignored any possible errors caused by the fitting of the stellar continuum. 

Errors in the derived temperatures and densities are dominated by two effects. Noise is 
an important source of error. For most PNe, we found that uncertainties due to this factor 
are about 5% in the derived temperatures and 0.1 dex in the densities. The other major 
error source originates from the extensive wings of the strong [O II] AA 3726, 3729 lines, 
which can affect the shape of the continuum near the Balmer jump. The effect is larger for 
low-excitation PNe, such as NGC 40 and NGC 6720. The extensive wings of the [O II] lines 
can cause the continuum redwards of the Balmer jump to be overestimated. As a result, 
the temperatures can be overestimated. The errors can be reduced by including data points 
at longer wavelengths in the fitting, where the extensive wings of the [O II] fines have less 
impact. However, increasing the wavelength range of the fitting would increase uncertainties 
in the reddening correction and in fitting the stellar continuum. To evaluate errors due to this 
effect, each spectrum was fitted twice for two wavelength ranges with one extending further 
into the red than the other. The extensive wings of the [O II] lines can also decrease the 
shape of the continuum-like feature redwards of the Balmer jump, leading to overestimated 
densities. For most PNe, this effect is negligible given the large wavelength difference (about 
70 A) between the Balmer jump and the moderate strength of [O II] fines. However, the 
effect can be significant for those with extremely strong [O II] fines. On the other hand, the 
extensive wings of the [O II] lines have much less effect on the integrated intensities of high- 
order Balmer lines {n 2, n = 11, 13, ...,23) which can be determined using a line profile 
fitting technique. Comparison of densities derived from the continuum-like feature redwards 
of the Balmer jump with those deduced from the integrated intensities of high-order Balmer 
lines (e.g. Liu et al. 2001b) allows us to estimate the errors in the derived densities caused 
by the extensive wings of the [O II] lines. The final errors given in Table 2 were based on 
estimates of contributions from the two dominating sources of uncertainty discussed above. 
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4. Discussion 

4.1. Comparison of densities 

Electron densities deduced from the Balmer decrements, and from the [O II] A3729/A3726, 
[S II] A6716/A6731, [CI III] A5517/A5537 and [Ar IV] A4711/A4740 doublet ratios are com- 
pared in Fig. 3. Fig. 3 shows that for most PNe, the hydrogen recombination spectrum 
yields higher densities compared to forbidden lines diagnostics. There is a trend such that 
the discrepancy becomes larger at higher densities. Two extreme cases are IC 4997 and 
M 2-24. For IC 3568, [CI III] A5517/A5537 yields an abnormally low density The cause 
of it is unclear but is probably due to measurement error. In Fig. 4 we compare electron 
densities derived from the [O III] far-infrared fine structure line ratio and from the hydrogen 
recombination spectrum. Evidently, densities derived from the [O III] 52/im/88/xm ratio 
are significantly lower than those derived from the Balmer decrements. In fact, it has been 
found that the far-IR hues yield systematic lower densities, lower than even those deduced 
from optical forbidden line density-diagnostics (Liu et al. 2001a). A linear fit shows that the 
density discrepancy increases with the electron density. 

0°, S°, C1+, 0+ and Ar^+ have ionization potentials of 13.6, 10.4, 23.8, 35.1 and 40.7 eV, 
respectively, and therefore embrace a wide range of regions of differing states of ionization. 
On the other hand, the hydrogen recombination spectrum arises from the entire ionized neb- 
ula. Thus the discrepancies between densities derived from H I recombination spectrum and 
those deduced from forbidden line ratios are not likely to be caused by ionization stratifica- 
tions. We notice that the [O II], [S II], [CI III], [Ar IV] and [O III] density-diagnostic lines 
have critical densities of iVcrit ^ 10^ cm^^. Therefore, if the nebulae contain regions of higher 
densities, these forbidden hues will be suppressed in these regions by coUisional de-excitation. 
In contrast, the hydrogen recombination spectrum is practically unaffected by coUisional de- 
excitation. As a result, in such cases the hydrogen recombination spectrum will yield higher 
average densities compared to the [O II], [S II], [CI III], [Ar IV] and [O III] lines ratios. The 
same argument has been developed by (Liu et al. 2001a) to explain the systematically lower 
densities derived from the [O III] fine-structure line ratio compare to those deduced from 
optical forbidden line ratios. The [Ar IV] lines have higher critical densities than the [O II], 
[S II], [CI III] lines. Thus, the effects are weaker for [Ar IV] doublet ratio than for the other 
cases. Given that the [O III] far-infrared fine-structure lines have the lowest critical densities 
amongst those density-diagnostic forbidden lines (iVcri([0 III]fr) ^ 4 x 10^ cm~^), the effects 
are largest for the [O III] 52//m/88//m ratio. This can be seen clearly in Fig. 3 and Fig. 4. A 
nebular model including these high-density condensations and their effects on temperature 
determinations will be discussed in Section 4.3. 



-14- 



4.2. Compcirison of temperatures 

Fig. 5 and Fig. 6 show that electron temperatures derived from the [O III] A4959/A4363 
and (52yum + 88;um)/A4959 ratios [hereafter, Tc([0 IIl]na) and Tc([0 IIl]fn)] are systematically 
higher than those deduced from the hydrogen recombination spectrum [hereafter Te(Bal)], 
although there are a few exceptions where Te(Bal) is higher than Te{[0 IIl]na). The most 
extreme case is He 2-118, for which T'e(Bal) is 5000 K higher than Te{[0 IIl]na)- In our 
sample, the four PNe with the largest temperatures discrepancies, Hf 2-2, M 1-42, M 2-36 
and NGC 6153, were previously studied in detail by Liu (2003), Liu et al. (2001b) and Liu 
et al. (2000). 

The discrepancy between Te([0 IIl]na) and re(Bal) was first discovered by Peimbert 
(1971) and was attributed to temperature fluctuations within the nebula (Peimbert 1967). 
He characterized the temperature structure of a nebula by the average temperature To and 
a mean square temperature fluctuation as follows: 

To AT^, Ni) = ' ' — 12 



and 



^2 _ /(Te - TofKNidV 



- TijN^NdV ' 

where Ni is the ion density for an observed emission line. Assuming that all the oxygen is 
twice ionized and that there are no density and composition variations in the nebula, Tq and 
can be derived from two of the three expressions, 

1 9 13x10^ 

Te([OIIlU=To[l + -(-^^ Z)t\ (14) 

re([o iii]f.) = ro[i + \{^^^^^ - (15) 

and 

re(Bal) =To(l-1.67i'), (16) 
(Peimbert 1967; Dinerstein et al. 1985). 

In Fig. 5 and Fig. 6, lines showing Tc([0 IIl]na) as a function of Tc(Bal) for the case of 
— 0.00, 0.02, 0.06 and 1.00 are over plotted. Fig. 5 shows that most nebulae have < 0.06. 
Excluding a few extreme nebulae having > 0.1, we obtain a mean value of < >= 0.031. 
The value is in excellent agreement with previous result < >— 0.03 derived by Liu & 
Danziger (1993) although their analysis contained larger measurement errors. However, 
even excluding the extreme nebulae, a similar comparison between T'e([0 Hljfn) and Te(Bal) 
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yields an average < t'^ >= 0.162. As we shall show in the following subsection, it is hkely that 
Te([0 Illjfn) have been overestimated because of density variations in the nebula, resulting 
in higher < t'^ > when comparing Te([0 Illjfn) and re(Bal). 

Another interesting point revealed by Fig. 5 is that Te(Bal) covers a wider range of 
values than Te([0 Illjna)- Except a few extreme cases, Te{[0 IIl]na) falls between 8000 and 
14000 K, whereas Te(Bal) varies from 900 to 19000 K. Peimbert (1967) has argued that tem- 
perature fluctuations of gaseous nebulae can cause lower Tc(Bal) compared to Te([0 IIl]na)- 
However, even < >= 0.031 is beyond any photoionization model prediction. Another 
explanation of higher Te(Bal) compared to Tc([0 IIl]na) is that there arc some condensa- 
tions with A'e > 10^ cm~^ in PNe, which then suppress the Tc([0 III]) AA4959,5007 relative 
to the A4363, leading to an artificially higher Te([0 Illjna) (Viegas & Clegg 1994). But no 
such high density condensations have been found in PNe having extremely large temperature 
discrepancies (Liu et al. 2000, 2001b). Recently, Liu et al. (2000) presented an empirical neb- 
ular model containing two components, each with its own temperature, density and chemical 
composition to account for the temperature discrepancy. This seems to be the most plausible 
explanation at present. 

For PNe having Te(Bal) > T'e([0 IIl]na), as in He 2-118, the electron temperature might 
conceivably increase in low-ionization regions owing to heating by shock waves in the outer 
regions. As a result, the average temperature of the entire nebula [Te(Bal)] appears higher 
than that of high-ionization regions [Te([0 Illjna,)]- 

4.3. A simple model with condensations 

In § 4.1, we found that electron densities derived from the H I recombination spectrum 

are systematically higher than those derived from optical forbidden lines ratios, and confirm 
the earlier result of Liu et al. (2001a) that densities derived from the [O III] far-IR fine- 
structure lines are the lowest of all. These results show that nebular electron densities 
derived from various diagnostics are correlated with the critical density of the diagnostic 
lines, and suggest that dense clumps are generally present in PNe. 

In such cases, the [O III] (52yum + 88;um)/A4959 ratio is no longer a good temperature 
indicator given the relatively low critical densities of the [O III] 52/im- and 88-/im lines 
(3500 and 1500 cm~^, respectively). When nebulae contain condensations with densities in 
the range between 10^ and 10^ cm~^, the [O III] 52/xm- and 88-/xm far-infrared lines will be 
heavily suppressed by coUisional de-excitation, leading to overestimates of Te{[0 Illjfn). On 
the other hand, condensations of densities in this range will not affect Te{[0 Illjna) given 
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that the [O III] A4959 and A4363 lines have higher critical densities of 7 x 10^ cm ^ and 
3 X 10^ cm~^, respectively (cf. Fig 6). 

Another piece of evidence pointing to the presence of dense clumps comes from Fig. 7 
where we compare T'e([0 IIl]na) and Te{[0 IIl]fn). If nebulae have temperature fluctuations 
but are homogeneous in density, then Te([0 III]fn) is expected to be always lower than 
Te([0 IIl]na) (Dinerstein et al. 1985). However, Fig. 7 shows about two-third of the objects 
have higher Tc([0 III]fn) compared to Tc([0 III]na), suggesting that density variations play 
an important role in Te([0 Illjfn) determinations. 

Viegas & Clegg (1994) showed that electron temperatures derived from the [O III] 
A4363/(AA4959 + 5007) ratio may be overestimated if nebulae contain ionized condensations 
of electron densities higher than 10® cm^'^. However, except for a few extreme PNe, such 
as IC 4997 (Hyung et al. 1994), Mz 3 (Zhang & Liu 2002) and M 2-24 (Zhang & Liu 
2003) there is no evidence that nebulae contain a substantial amount of ionized gas in 
condensations with densities in excess of 10^ cm~^, as the presence of such condensations 
should have been revealed by the H I recombination spectrum (Liu et al. 2000, 2001a, 2001b; 
Fig. 3). Therefore, we consider A^e ~ 10^ cm~^ as a safe upper hmit for the density of these 
condensations. 

In order to estimate the filling factor of the condensations, we use an analytical method 
originally developed by Viegas & Clegg (1994) for explaining the discrepancy of T'e([0 IIl]na) 
and re(Bal). Assuming that some high-density regions {N^ > 10^ cm~^) are present, they 
deduced nebular-to-auroral [O III] line ratios (see their Eq. 2.3). Following the same deduc- 
tion of their Eq. 2.3 but assuming that the electron density of these condensations is about 
10^ cm~^ instead of 10® cm~^, we obtain the infrarcd-to-auroral line ratio 

-^52/Xm+88//m _ £fL 1 + /i^Cj(££H/gfL) ^^^^ 
-^4959 £nL 1 + Af^<^ (Snu/ ^riL ) ' 

where e is line emission coefficient, AttJ = ncUie, f and n refer to IR fine-structure line 
and nebular line, L and H refer to the low- and high-density regions, /x = (A''c)H/(A''e)L is 
the density contrast and uj = Vh/Vl is the filling factor. In the following calculation, we 
assume that: 1) density variations are the only cause of the lower Te([0 Illjfn) compared to 
T'e([0 IIl]na) and Te(Bal); 2) the nebular gas has a homogeneous electron temperature, i.e., 
re([0 IIl]na) Or Te(Bal); 3) the condensations have densities lower than the critical density 
of the [O III] auroral fine but higher than the critical densities of the [O III] infrared lines, 
whereas low-density regions have a density that is lower than the critical densities of the 
[O III] auroral and infrared lines. 

Setting ^(^[o ni]f„) = -^52//m+88//m/-^4959, we have £fL/£nL = -^(Tsai)- For a simple 
assumption (Ae)H = 10^ cm~^ and (Ae)L = 10^ cm~^, the ratio Sfn/^fL is about 0.04, 



whereas £„ is practically independent of the density when the densities are lower than the 
critical density of the nebular line and thus Snn/SnL ~ 1- Therefore, /i^cu is given by 




F(T[o in],J/F(TBai) - 1 



(18) 



0.04-F(T[o mi,J/F(TBai)- 



In Fig. 6, the hues showing the variation of T'e([0 IIl]fn) as a function of Te(Bal) for 
fx^uj — 0.0, 0.1, 0.5, 1.0 and 2.0 are plotted. We obtain an average value of < /ji^u; >— 1.0. 
However, it should be borne in mind that the n'^uj derived by eq. (18) is a upper limit 
since temperature variations are actually present and contribute partly to the measured 
temperature discrepancies. 

fi'^uj can also be estimated by comparing Te([0 Illjna) with Te([0 IIl]fn), as shown in 
Fig. 7. Assuming that T'e([0 Illjna) is constant across the entire nebula, /s^cu is likewise given 



In Fig. 7, the lines showing the variation of Te([0 Illjfn) as a function of Te(Bal) for ii'^uj = 0.0, 
0.1, 0.5, 1.0 and 2.0 are also plotted. An average value of < /I'^u; >— 0.4 is obtained. It 
is noteworthy that the presence of temperature variations will increase Te([0 Illjna) with 
respect to T'e([0 Illjfn), in contrast to density variations. Therefore, //^o; derived from eq. (19) 
is a lower limit. 

The upper- and lower-hmit of /jl^uj are hsted in Table. 3. For a few PNe having a lower 
limit (/i^a;)L = 0.0, the effect of temperature variations on the temperature determination is 
larger than that of density variations, causing higher Te{[0 Illjna) compared to Te([0 Illjfn). 
For all the PNe, the mean value of /I'^cu is in the range from 0.4 to 1.0. Thus, we suggest that 
n'^LV = 0.7 may be a representative value that can be used for a typical nebula. Adopting 
the density contrast fi ~ 100, a representative value of the condensations filling factor uj is 
about 7 X 10^^. The value is similar to that found by Viegas & Clegg (1994). Adopting 
/i^a; = 0.7 for all the PNe, we obtain new Te([0 Illjfn), for which the contribution of the 
clumps has been deduced using equation (17). Comparing the new Te([0 Illjfn) and Te(Bal), 
we re-evaluate the average value < t'^ >— 0.032, which is in good agreement with the value 
yielded by the comparison of T'e([0 Illjna) and re(Bal) (see Section 4.2). 



by 




FjTlO III]f J/-^(^[0 III]na) ~ 1 



(19) 



0.04-F(T[o iii],J/F(T[o iii]„J' 
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4.4. The difference between temperatures derived from the [O III] forbidden 
Unes and from the Balmer jump versus electron density 

In Fig. 8, the differences between temperatures derived from tlie [O III] forbidden lines 
and from the Balmer jump are plotted against electron densities derived from the hydrogen 
recombination spectrum, the [O II] A3729/A3726, [S II] A6716/A6731, [CI III] A5517/A5537 
[Ar IV] A4711/A4740 and the [O III] 52//m/88//m ratios, respectively. In all cases, the 
differences between the two temperature are anti-correlated with nebular densities derived 
from various density indicators. A least-squares fit yields the following relation: 

AT= (1.30 ±0.12) X 10^ - (2.93 ±0.27) x lO^logiVe, (20) 

and a hnear correlation coefficient of 0.51, where AT = Te([0 III])— Te(Bal) and A^e is that 
derived from hydrogen recombination spectrum, which represents the average density of 
the entire nebula. Garnett & Dinerstein (2001) found that the discrepancies of O^^/H"*" 
abundance ratios derived from ORLs and from CELs are greater for larger, lower-surface 
brightness PNe than for compact, dense ones. On the other hand, Liu et al. (2001b) found 
that the abundance discrepancies are greater for those with larger temperature discrepancies. 
Our results are therefore consistent with these earher findings. Liu et al. (2001b) found that, 

A(02+/H+) = (0.209 ± 0.085) ± (2.01 ± 0.33) x 10"^ AT, (21) 

where A(0^+/H+) = log(0^+/H+)oRL — log(0^+/H+)cEL- From the two relations given in 
equations (20) and (21), we find that the ratio of abundances derived from ORLs and from 
CELs has a power-law dependence on density, 

q2+ 

^ORL ^ jy-a^ (22) 



q2+ ^^'e 



where a = 0.59 ± 0.15. 



In general, electron densities of PNe are expected to decline monotonically as nebulae 
expand. Schmidt- Voigt & Koppen (1987) found an approximate relation between nebular 
mean electron density and nebular size, A'e oc R~^, where R is the nebular radius. Conse- 
quently, it is suggested that the temperature discrepancy is related to the evolutionary state 
of PNe; more expended PNe show larger temperature discrepancies. On the other hand, 
the actual evolution in a given nebula may depend on many factors. Furthermore, the elec- 
tron density measured for a given nebula may depend on the diagnostic used. For example, 
electron densities of low- ionization regions, as probed by the [O II] A3729/A3726 and [S II] 
A6716/A6731 line ratios, can be enhanced by interactions of fast stellar winds with earlier 
slow winds. Therefore, there are shght difference in deriving the relation of AT and from 
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different density diagnostics, as shown in Fig. 8. We also use electron density derived from 
these forbidden line ratios to estimate the a- value in equation (22). The results show that 
the value is within the range of 0.3 and 0.9. 

To explain the relation shown in equation (20) and (22), we apply a chemically inho- 
mogeneous two-component model presented by Liu et al. (2000). In their hypothesis, the 
nebula consists of two components, H-deficient clumps and diffuse material with 'normal' 
abundances (~solar). These H-deficient clumps are expected to have a low electron temper- 
ature due to enhanced cooling caused by high metallicity and pressure equilibrium between 
the two components, so that only recombination lines arc emitted. As a result, ORLs yield 
overestimated heavy element abundances. In a given PN, the contribution of the two compo- 
nents to the Balmer temperature and abundance determinations depends on their densities, 
i.e., the contribution becomes weaker as density declines. If the clumps have neutral cores, 
as found in Abell 30 (Ercolano et al. 2003) and NGC 7293 (Dyson et al. 1989), they are 
expected to be quite stable during the evolution of the PNe. As a result, the relative con- 
tribution of the low-density and high-temperature material to the average temperature and 
abundance determinations decrease as nebulae expand. Consequently, the discrepancies in 
temperature and abundance determinations become larger as the nebulae expend. 



4.5. Are He/H ratios generally overestimated? 

A long-standing problem in abundance determinations for PNe has been that the heavy 
clement abundances derived from ORLs are systematically higher than those derived from 
CELs. The ratio of abundances derived from the two types of emission lines is found to vary 
from target to target. Tsamis et al. (2003a) analysed a PN sample and found the ORL/CEL 
discrepancy factors for these heavy element abundances span a range from about 2 to as large 
as 20. For the most extreme known PNe Hf 2-2, the ORL/CEL abundance ratio reaches a 
record value of 84 (Liu 2003). Another problem in nebular astrophysics is that T'e([0 IIl]na) 
are systematically higher than T'e(Bal), as we have discussed in § 4.2. The two problems are 
found to be related. A positive correlation between the ratio of 0^'^/H"'" derived from O II 
ORLs and from [O III] CELs and the temperature discrepancy [Ar=Te([0 IIl]na)-re(Bal)] 
has been shown by Liu et al. (2001b). 

The explanations for the dichotomy in abundance and temperature determinations have 
focussed on temperature and density variations and chemical inhomogeneities. However, it 
has been shown that the temperature and density variations alone can be ruled out by current 
observations (Liu 2001, 2003). At present, the most plausible explanation is that there exist 
some cold H-deficient clumps in nebula (Liu et al. 2000) where only recombination lines are 
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emitted. 

To fit the observed pattern of NGC 6153, Liu et al. (2000) liave presented several em- 
pirical models incorporating H-deficient clumps. Pequignot et al. (2003) also constructed 
photoionization models of NGC 6153 and M 1-42, incorporating H-deficient clumps , which 
reproduce the observed integrated line ratio satisfactorily. In all the models, the H-deficient 
clumps contain only about 1% of the total mass, suggesting that the H-deficient components 
contain only a tiny amount of nebular gas. Nevertheless, they can seriously enhance the in- 
tensities of ORLs, causing overestimated ORL abundances. Therefore, the He/H abundance 
ratio could be overestimated owing to the existence of the H-deficient clumps can lead to 
since all the observable strong He I and He II lines are recombination lines. Some H-deficient 
components have been found in some PNe, such as Abell 30 and Abell 78 (Jacoby & Ford 
1983). Wesson et al. (2003) analysed optical spectra of the H-deficient knots in Abell 30 and 
found the He/H abundance ratios of the knots are more than ten. 

These PNe with extremely large abundance and temperature discrepancies, such as 
NGC 6153, M 1-42 and Hf 2-2, have been found to have high He/H abundances. Peimbert & 
Torres-Peimbert (1983) defined the PNe having high He/H abundance (> 0.125) as Type-I 
PNe. According to our suggestion, these PNe with high He/H abundances might result 

from contamination from H-dcficient clumps. Quantitative analysis of overestimated He/H 
abundances in PNe is quite difficult. For this purpose, as mentioned by Pequignot ct al. 
(2003), one has to rely on the measurement of He I A10830 line which is strongly enhanced 
by collisional excitation from the 2s^S metastable level. 

If the He/H abundance is enhanced by H-dcficient clumps, one expects that there is 
a positive correlation between the He/H abundance and the ORL/CEL abundance ratio of 
heavy element. For most PNe, however, heavy element ORLs are quite weak and thus yield 
large uncertainties in the resultant abundance. Instead, we can study the relation between 
the He/H abundance and AT since the ORL/CEL abundance ratio has been found to be 
positively correlated with AT by Liu et al. (2001b). Fig. 9 plots He/H versus AT. The 
measurement error of the He/H abundance ratio is negligible. In the figure, filled and open 
circles represent the Galactic disc and bulge PNe respectively and the solid line is a hnear 
fit represented by the relation: 

He/H = (0.056 ± 0.003) + (4.25 ± 0.11) x 10"^ AT, (23) 

which yields a linear correlation coefficient of 0.50. Although there is considerable scatter, 
there is a trend of increasing He/H abundance with increasing difference between the tem- 
perature from [O III] forbidden fines and Balmer discontinuity for both the Galactic disc 
and bulge PNe. Therefore, it seems plausible that the He/H abundances in PNe have been 
generally overestimated due to the existence of cold H-deficient clumps. 
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Fig. 9 also shows that the most of PNe have enhancement of hehum with respect to 
the Sun, which is often ascribed to the second and third dredgc-up. However, the normal 
model of dredge- up cannot explain some extreme cases, such as He 2-111 which has a He/H 
abundance of 0.219 (Kingsburgh & Barlow 1994). For these cases, the more likely reason 
is that pollution of the H-deficient clumps causes an overestimate of the He/H abundance 
ratio. On the other hand, we cannot completely reject the possibility that helium in nebulae 
has been partly enhanced by the products of nucleosynthesis of the central star. That may 
be the main cause of the scatter in the derived He/H abundances. The presence of hydrogen- 
deficient clumps may be further proof. 

The ORL/CEL abundance discrepancies are also found in H II regions (Esteban ct al. 
2002; Tsamis et al. 2003b). And the existence of some high density clumps in H II regions 
has been indicated by previous studies (e.g. O'Dell et al. 2003). Therefore, the H-deficient 
material might also be present in H II regions, leading to overestimated ORL abundances. 
If this is the case, the primordial helium abundance, which is based on observations of some 
metal poor extragalactic H II regions, could have been overestimated. 

5. Temperatures and densities derived from the hydrogen recombination 

spectra near the Paschen jump 

The hydrogen recombination spectrum near the Paschen jump can also provides a di- 
agnostic of electron temperature and density [hereafter, Te(Pas) and A^'e(Pas)]. Unfortu- 
nately, in our sample, the spectra covering the Paschen jump are available for only four 
PNe, NGC 7027, NGC 6153, M 1-42 and NGC 7009. Using the same method described in 
Section 3, wc match the theoretical with the observed spectra at 8250 A. Given that the 
observed spectra near the Paschen discontinuity may include a contribution from a dust 
emission tail, we added a constant continuum for the match. The results are shown in 
Fig. 10. All the spectra are normahzed to P 20 in order to decrease the uncertainty of red- 
dening correction. Table 4 lists the derived temperatures and densities. One advantage of 
using the spectrum near the Paschen jump is that the resultant temperature and density are 
not affected by the continuum of the hot central star and two-photon emission. However, 
the emission from local dust may be stronger in this band. In addition, the estimation of 
the continuum blucwards the Paschen jump is affected by the telluric absorption (H2O and 
O2). And the Ca II absorption feature might also cause a error in fitting the continuum 
longward of the Paschen jump of some PNe, such as M 1-42. As a result, the temperatures 
and densities derived from the hydrogen recombination spectra at the Paschen jump are less 
reliable than those derived from those at the B aimer jump. 
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Table 4 shows that the Ne{Pets) are consistent with the A^e(Bal) within the errors. 

However, the Te(Pas) are found to be 4000, 1000, 500 and 1200 K lower than the Te(Bal) for 
NGC 7027, NGC 6153, M 1-42 and NGC 7009, respectively. This may be due to systematic 
errors in the reddening correction or flux cahbration. As an ahernative, it remains possible 
that the discrepancy is caused by temperature variations. Further investigations including 
a larger data sample of PNe is needed to clarify the point. 



6. Summary 

In this work, we presented a method to derive nebular electron temperature and density 
simultaneously by matching the theoretical with the observed hydrogen recombination spec- 
tra. Using the method, we derived electron temperatures and densities of 48 Galactic PNe. 
The temperatures and densities derived by this method are compared with those derived 
from other indicators, enabling us to study electron temperature and density variations in 
PNe. Our main findings are as follows. 

— The densities derived from hydrogen recombination spectra are generally higher than 
those deduced from forbidden line ratios, suggesting that condensations are generally 
present in nebulae. 

— The temperatures derived from the nebular- auroral line ratios of [O III] are system- 
atically higher than those derived from hydrogen recombination spectra, suggesting 
that temperature variations are generally present in nebulae. We obtained the rms 
temperature fluctuation parameter t'^ — 0.031 as a representative value. 

— The larger discrepancy in temperatures derived from [O III] (52/im -|- 88/im)/A4959 
line ratios and from hydrogen recombination spectra, however, cannot be attributed 
simply to the temperature variations. The presence of clumps with A^e ~ 10^ cm~^ 
can account well for the discrepancy. A comparison between electron temperatures 
derived from [O III] (52/im -|-88/im)/A4959 fine ratios and those derived from hydrogen 
recombination spectra and [O III] A4959/A4363 line ratios suggests that the flUing 
factor of the condensations has a typical value of 7 x 10~^. 

— The discrepancies in temperatures derived from [O III] A4959/A4363 hne ratios and 
from hydrogen recombination spectra are related to the evolution of PNe. The ORL / GEL 
abundance ratio is found to have a power-law dependence on electron density (~ A^", 
a = 0.59±0.15). We attribute this to different evolutionary scenario of two components 
in PNe, i.e. expanding diffuse nebular material and stable cold H-deficient clumps. The 
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contribution of the clumps to the temperature and abundance determinations increases 
as the nebulae expand. 

— He/H abundance has a positive correlation with the temperature discrepancy, suggest- 
ing that He/H abundances might be overestimated due to the existence of H-deficient 
clumps. If this is the case, traditional up-dredge theory should be modified. The 
primordial helium abundance determination will be affected if there are significant 
condensations in H II regions. Further investigations are needed. 

— The densities derived from the spectrum near the Paschen jump are in good agree- 
ment with those derived form that near the Balmer jump. However, the former yield 
generally lower temperatures than the latter. The reason remain unknown. 

We conclude that the existence of H-deficient clumps in PNe is the most plausible ex- 
planation for the origin of temperature and density variations and the ORL/CEL abundance 
discrepancies. However, it is a challenge to clarify the origin and evolution of these clumps. 
In further work, wc will study the statistical relationships between the temperature and 
abundance discrepancies and other physical conditions of PNe such as radio flux, morphol- 
ogy, chemical abundances and the nature of the central star. 
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Fig. 1. — Theoretical hydrogen recombination spectra near the Balmer discontinuity at (a) 
same densities but different temperatures and (b) same temperatures but different densities. 
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Fig. 2. — A comparison of observed and theoretical spectra at 3650 A. The top of each panel 
shows the observed spectrum (solid line) which has been de-redddened using the extinction 
constant listed in Table 1, the synthetic spectrum (dotted line), nebular continuum (dotted- 
dashed line) and stellar continuum (dashed hne) . The bottom shows residuals subtracted by 
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log Ne(Bal) (om-3) log Ne(Bal) (om-3) 

Fig. 3. — Comparison of the electron densities derived from the intensity ratios of [O II] 
A3729/A3726, [S II] A6716/A6731, [CI III] A5517/A5537 and [Ar IV] A4711/A4740 versus 
those deduced from the Balmer decrements. The sohd hues are y — x plots. 
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Fig. 4. — Comparison of the electron densities derived from [O III] 52/im/88/i.m ratios and 
the Balmer decrements. The sohd hne is a y — x plot and the dashed line is a linear fit. 
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5000 10000 15000 20000 

Te(Bal) (K) 

Fig. 5. — Comparison of the electron temperatures derived from the nebular-to-auroral 
line ratio of [O III] and those deduced from the Balmer decrements. The lines show the 
variation of Te{[0 IIl]na) ^ ^ function of Te(Bal) for the mean square temperature fluctuation 
= 0.00, 0.02,0.06 and 0.10. 
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Fig. 6. — Comparison of the electron temperatures derived from the [O III] (52/im + 
88/im)/A4959 ratios and those deduced from the Balmer decrements. The dotted hues show 
variation of Te([0 Illjfn) as a function of Te(Bal) for the mean square temperature fluctua- 
tion t'^ — 0.02, 0.06 and 0.10. The dashed lines show variation of Te([0 Illjfn) as a function 
of re(Bal) for pi'^u; = 0.1,0.5,1.0 and 2.0 (see Section 4.2). The solid line shows that for 
homogeneous temperature and density. 
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Fig. 7. — Comparison of the electron temperatures derived from the [O III] A4959/A4363 
ratios and those deduced from the [O III] (52yum+88/im) / A4959 ratios. The dotted hues show 
the variation of Te([0 III]na) as a function of Te([0 III]fn) for the mean square temperature 
fiuctuation — 0.02,0.06 and 0.10. The dashed hnes show the variation of Te([0 Illjna) 
as a function of T'e([0 Illjfn) for fi'^u; — 0.1,0.5, 1.0 and 2.0. The solid line shows that for 
homogeneous temperature and density. 
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Fig. 8.- Te([0 IIl])-Te(Bal) versus 7Ve(Bal), N,{[0 II]), N,{[S II]), 7V^([C1 III]), 7V^([Ar IV]) 
and A^e([0 IUff])- Tlie solid lines are least-squares fits. The linear correlation coefficient "r" 
is given in every panel. For the fits the observed points are vs^eighted according to the error 
bar for each. For all the trend of increasing discrepancy of temperature derived from 

[O III] forbidden lines and the Balmer jump with decreasing electron density is seen. 
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Te([0 III]) - Te(Bal) 

Fig. 9. — He/H versus Te([0 III])— TelBal). The filled and open circles are the Galactic disc 
and bulge PNe, respectively. The solid line is a least-squares fit. The dashed line is the solar 
He/H abundance. A trend of increasing He/H abundance with increasing difi^erence between 
the temperature from [O III] forbidden lines and Balmer discontinuity is seen, indicating 
that He/H ratios in PNe may be overestimated generally. 
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Fig. 10. — A comparison of 4 observed and theoretical spectra at 8250 A. The top of each 
panel shows the observed spectrum (solid line) which has been deredden using the extinction 

constant listed in Table 1 and the synthetic spectrum (dot line) . The bottom shows residuals 
subtracted by 0.05. Note the absorptions near 8200 A by atmospheric gases (H2O and O2) 
and two Ca II absorption features at 8500 A and 8546 A in the spectrum of M 1-42. 
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Table 1. Object and parameters; line fluxes in units of 10 erg cm ^ s ^. 



Source 


c He+/H+ He2+/H+ Te([0 IIl]na) 


Te([0 IIl]fn) (To 




(K) 


(K) (Km/s) 



Hu 1-2 


0. 


,512 


0, 


,046 


0, 


.079 


19500 ± 300 




144 


TC 3568 


n 


,261 


n 


095 


n 


001 

• \J\J ± 


1 1 370 -1- 200 


12600 ± 200 


132 


NGC 7027 


1. 


380 


n 


058 


n 


041 


12600 + 180 

X^VJWW —1— xow 


1 5800 + 3000 


82 


NGC 700Q 





130 





093 





01 1 




8350 + 1 50 


132 


NGC 40 


0. 


,700 


0. 


,062 


0, 


.000 


10600 ± 200 


10000 ± 200 


132 


NGC 6210 





055 





099 





001 




1 nmc) -\- 200 


148 


NGC 6826 


n 


064 


n 

yj . 


099 


n 


000 


9350 -1- 200 


1 0600 + 200 


132 


NGC 6572 


0. 


,400 


0, 


,110 


0, 


.000 


10600 ± 200 


17100 ± 200 


132 


NGC 6741 


1, 


,100 


0. 


,081 


0, 


.031 


12600 ± 200 


19000 ± 2000 


132 


NGC 6790 


1, 


100 

■ ±\J\J 


n 


081 


n 

yj . 


003 


12840 + 200 




132 


NGC 6884 


1. 


050 


n 


081 


n 


016 






132 


NGC 7662 





180 





061 

. \J\J X 





038 


13300 ± 200 


13000 ± 200 


156 

xu w 


NGC 6543 





100 


n 


117 


n 


000 


7040 + 200 


8520 + 1 50 


140 


M 2-24 


n 


800 





099 


n 


005 


01 nno -1- 1 500 




156 


L>n z-i 


1 

i. 


,UDo 


U. 


1 H/l 
, iU4 


n 

U, 


.UU4 


iU4UU ± OUU 






NGC 6720 


0. 


,200 


0. 


,094 


0, 


.017 


10630 ± 200 


11600 ±200 


132 


NGC 6620 


0. 


,760 


0. 


,099 


0, 


.021 


10400 ± 400 




132 


NGC 6567 


0. 


,670 


0. 


,105 


0, 


.007 


11400 ±400 


13800 ± 300 


132 


H 1-35 


1. 


,360 


0, 


,115 


0, 


.001 


10500 ±500 




132 


H 1-50 


0. 


,700 


0. 


,100 


0, 


.011 


11300 ±500 




132 


He 2-118 


0. 


,158 


0. 


,065 


0, 


.000 


13400 ± 300 




148 


M 1-20 


1. 


,200 


0. 


,094 


0, 


.000 


10700 ± 500 




123 


M 3-21 


0. 


,570 


0. 


,100 


0, 


.007 


11100 ± 100 




140 


M 3-32 


0, 


,579 


0, 


,120 


0, 


.013 


9200 ± 200 




140 


IC 1297 


0. 


,220 


0. 


,079 


0, 


.042 


10700 ± 100 




132 


IC 4634 


0, 


,454 


0. 


,097 


0, 


.000 


10200 ± 100 


12500 ± 300 


140 


IC 4776 


0, 


,265 


0. 


,088 


0, 


.000 


11000 ±200 




136 


IC 4997 


0. 


,438 


0. 


,139 


0, 


.000 


20000 ± 1500 




127 


NGC 5873 


0. 


,004 


0. 


,063 


0, 


.040 


12800 ± 100 




127 
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Table 1 — Continued 



Source 




c 


rle / rl^ 


tie ^ / ti^ 


T" fir^ tttI ^ 

le([(J Illjnaj 


T" fir^ tttI ^ 

le([U illjfnj 
















(K) 


[ixm/s) 


NGC 5882 


0. 


,380 


0.113 


0.003 






127 


NGC 6302 


1. 


,460 


0.105 


0.067 


18900 ± 200 


15500 ± 500 


132 


NGC 6818 


0. 


,400 


0.050 


0.051 


12600 ± 200 


15700 ± 300 


148 


NGC 6153 


1. 


,300 


0.123 


0.011 


9110 ± 200 


8380 ± 200 


148 


M 2-36 


0. 


,270 


0.133 


0.003 


8380 ± 200 


8380 ± 160 


148 


M 1-42 


0. 


,700 


0.139 


0.009 


9220 ± 300 


6870 ± 200 


140 


NGC 6778 


0, 


,890 


0.155 


0.008 


9200 ± 300 




140 


Hf 2-2 


0. 


,200 


0.131 


0.002 


8820 ± 500 




132 


DdDm 1 


0. 


,137 


0.060 


0.000 


12330 ± 700 




99 


Hu 2-1 


0. 


,777 


0.097 


0.000 


9850 ± 300 




99 


IC 2003 


0. 


,347 


0.054 


0.049 


12640 ± 400 




99 


IC 5217 


0. 


,501 


0.088 


0.008 


11260 ±400 




115 


Me 2-2 


0. 


,343 


0.138 


0.000 


10960 ± 300 




90 


NGC 6803 


0. 


,869 


0.112 


0.004 


9730 ± 300 




107 


NGC 6807 


0, 


,642 


0.084 


0.000 


10920 ± 500 




107 


NGC 6833 


0. 


,000 


0.078 


0.000 


12810 ±400 




99 


NGC 6879 


0. 


,401 


0.091 


0.003 


10400 ± 600 




99 


NGC 6891 


0. 


,287 


0.094 


0.000 


9320 ± 300 




115 


Vy 1-1 


0. 


,421 


0.094 


0.000 


10240 ± 600 




99 
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Table 2. Electron temperatures and densities derived from the spectra near the Balmer 

jump and forbidden hne ratios. 



Source 


Te(Bal) 


logNc(Bal) 




(K) 


fcm~^) 


Hu 1-2 


20000 ± 1200 


3.9 ±0.3 


IC 3568 


9300 ± 900 


3.8 ±0.2 


NGC 7027 


12000 ± 400 


5.1 ±0.1 


NGC 7009 


7200 ± 400 


3.8 ±0.1 


NGC 40 


7000 ± 700 


3.2 ±0.3 


NGC 6210 


8700 ± 800 


3.8 ±0.1 


NGC 6826 


8700 ± 700 


3.4 ±0.1 


NGC 6572 


10300 ± 1000 


4.4 ± 0.2 


NGC 6741 


15000 ± 600 


3.6 ±0.2 


NGC 6790 


14000 ± 500 


4.5 ±0.1 


NGC 6884 


11600 ±400 


3.8 ±0.1 


NGC 7662 


12200 ± 600 


3.4 ±0.3 


NGC 6543 


6800 ± 400 


3.8 ±0.2 


M 2-24 


16000 ± 2000 


7.0 ±0.7 


Cn 2-1 


10800 ± 600 


4.8 ±0.2 


NGC 6720 


8000 ± 600 


2.9 ±0.3 


NGC 6620 


10000 ± 500 


3.7 ±0.2 


NGC 6567 


12000 ± 800 


4.0 ±0.5 


H 1-35 


10000 ± 1000 


4.8 ±0.2 


H 1-50 


12500 ± 1000 


4.1 ±0.3 


He 2-118 


18300 ± 1000 


5.0 ±0.5 


M 1-20 


12000 ± 700 


4.2 ±0.3 


M 3-21 


11000 ±400 


4.2 ±0.1 


M 3-32 


4500 ± 500 


3.6 ±0.2 


IC 1297 


10000 ± 400 


4.0 ±0.2 


IC 4634 


8500 ± 400 


4.3 ±0.2 


IC 4776 


8600 ± 400 


4.8 ±0.2 


IC 4997 


10200 ± 500 


6.0 ±0.2 


NGC 5873 


12000 ± 500 


3.9 ±0.1 


NGC 5882 


6800 ± 500 


4.0 ±0.2 
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Table 2 — Continued 



Source 




logNe^hsalj 






/',-„-, -,-3 A 

(cm j 


NGC 6302 


16100 ± 500 


4.9 ± 0.1 


NGC 6818 


12500 ± 600 


3.7 ± 0.1 


NGC 6153 


6000 ± 400 


3.8 ±0.2 


M 2-36 


6000 ± 400 


3.8 ±0.1 


M 1-42 


4000 ± 600 


3.7 ±0.2 


NGC 6778 


4300 ± 500 


3.6 ±0.2 


Hf 2-2 


1000 ± 400 


2.6 ±0.5 


DdDm 1 


11400 ±800 


3.8 ±0.2 


Hu 2-1 


10000 ± 500 


4.3 ±0.2 


IC 2003 


11000 ± 800 


3.0 ±0.5 


IC 5217 


12000 ± 800 


3.0 ±0.5 


Me 2-2 


11000 ± 500 


4.2 ± 0.2 


NGC 6803 


8500 ± 400 


3.6 ±0.4 


NGC 6807 


12000 ± 1000 


3.5 ±0.6 


NGC 6833 


14000 ± 500 


4.1 ±0.3 


NGC 6879 


7000 ± 400 


3.3 ±0.3 


NGC 6891 


6000 ± 600 


3.6 ±0.6 


Vy 1-1 


6000 ± 600 


3.0 ±0.5 



-40- 



Table 3. The evaluation of ji'^u. 



Source 




F(T\o III],) 


F{T\o iiii„„) 


(//^a;)H 




IC 3568 


0.527 


0.254 


0.315 


1.174 


0.253 


NGC 7027 


0.281 


0.169 


0.254 


0.708 


0.532 


NGC 7009 


1.210 


0.720 


0.433 


0.730 


0. 


NGC 40 


1.345 


0.433 


0.374 


2.408 


0. 


NGC 6210 


0.641 


0.392 


0.469 


0.681 


0.208 


NGC 6826 


0.641 


0.374 


0.520 


0.768 


0.415 


NGC 6572 


0.401 


0.151 


0.374 


1.857 


1.639 


NGC 7662 


0.271 


0.238 


0.228 


0.145 


0. 


NGC 6543 


1.505 


0.688 


0.872 


1.300 


0.282 


NGC 6720 


0.836 


0.303 


0.374 


1.981 


0.246 


IC 4634 


0.688 


0.258 


0.411 


1.867 


0.634 


NGC 6818 


0.258 


0.171 


0.254 


0.543 


0.516 


NGC 6153 


2.543 


0.714 


0.254 


2.987 


0. 


M 2-36 


2.543 


2.056 


0.714 


0.249 


0. 


M 1-42 


2.543 


1.422 


0.543 


0.850 


0. 


Average 








1.0 


0.4 
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Table 4. Electron temperatures and densities derived from the spectra near the Paschen 

jump. 



Source 


Te(Pas) 


logNe(Pas) 




(K) 


(cm~^) 


NGC 7027 


8000 ± 1000 


5.4 ±0.7 


NGC 6153 


5000 ± 700 


3.8 d= 0.5 


M 1-42 


3500 ± 500 


3.0 ±0.8 


NGC 7009 


5800 ± 900 


3.5 ±0.7 



